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Abstract 

We put forward the concept of quantum spiral bandwidth of the spatial mode function of the 
two-photon entangled state generated in spontaneous parametric down-conversion. We obtain 
the bandwidth using the eigenstates of the orbital angular momentum of the biphoton states, and 
reveal its dependence with the length of the down-converting crystals and waist of the pump beam. 
The connection between the quantum spiral bandwidth and the entropy of entanglement of the 
quantum state is discussed. 

PACS numbers: 03.67.-a, 42.50.Dv 
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Entanglement, one of the most genuine features of quantum mechanics, is a basic ingredi- 
ent in quantum cryptography, computing and teleportation □ D . Spontaneous parametric 
down-conversion, namely the generation of two lower frequency photons when a strong 
pump field interacts with a nonlinear crystal, is a reliable source of entangled photons. The 
generated two-photon state is entangled in polarization [3], and most applications ofpara- 
metric down-conversion in quantum systems make use of such spin-entanglement Ql-j^l- 
However, entanglement is also embedded in the spatial mode functions that describe the 
two-photon states Such spatial entanglement occurs in an infinite dimensional Hilbert 
space and is gaining increasing attention as a powerful way to encode and to exploit quantum 
information. |a S llfl llll, Q, 13 For example, knowledge of the spatial structure acquired by 
the mode functions of entangled signal and idler photons forms the basis of multidimensional 
quantum imaging, and it can be used to increase the efficiency of multidimensional quan- 
tum communication protocols. The spatial structure of the photon states can be expressed 
by a mode decomposition of their mode function in an appropriate basis. The amount of 
entanglement of a quantum state is directly related to the width of such modal expansion, 
hereafter referred to as quantum spatial bandwidth. A two-photon state described by a single 
mode is a product state, while a quantum state described by an equi-distributed multimode 
expansion corresponds to a maximally entangled state. Therefore, a fundamental question 
that arises is how to act on the spatial quantum distribution of a given mode expansion, 
hence how the corresponding quantum information can be, e.g., concentrated. 

It has been recently observed experimentally that the photon pairs generated in sponta- 
neous down-conversion are entangled in orbital angular momentum (OAM) Such OAM 
entanglement allows the implementation of arbitrary d- dimensional quantum channels jl^ . 
thus it has been used to demonstrate violation of Bell inequalities with qutrits The 
corresponding mode functions are naturally expressed in terms of eigenstates of the paraxial 
OAM operator, whose spiral or winding topological structure can be resolved experimentally 
using combinations of holographic and filtering techniques. Therefore, such base provides a 
powerful tool to explore the concept of spatial multimode entanglement. In this Letter we 
put forward the full characterization of the entangled photon pairs in terms of eigenstates of 
the OAM operator and reveal how the corresponding quantum spiral bandwidth depends on 
the shape of the beam that pumps the down-converting crystal, and on the material proper- 
ties and length on the crystal. We obtain in analytical form important similarity rules that 
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hold for arbitrary input and crystal conditions. 

Let a quadratic nonlinear crystal of length L be illuminated by a quasi-monochromatic 
laser pump beam propagating in the z direction. The spatial distribution of the pump beam 
writes E p (x, z, t) = E (x) exp [i {uj p t — k p z)} + c.c, where u p is the angular frequency of the 
pump beam, k p is the wave number, x is the position in the transverse plane and E is the 
transverse spatial field profile at the input face of the crystal. The length of the crystal is 
assumed to be much smaller than the diffraction length (Ld) of the pump beam (L << La), 
so that the spatial distribution of the pump beam in the transverse plane is independent 
of z. The signal and idler photons are assumed to be monochromatic, with uj p = uo s + Ui, 
where and uj s and are the frequency of the signal and idler photons, respectively. This is 
justified by the use of narrow band interference filters in front of the detectors. 

In the paraxial approximation, the spin and the OAM can be considered separately jl7|. 
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Under these conditions, photons described by a mode function that is a 
(LGp) mode are eigenstates of the OAM operator with eigenvalue Xh 
is the number of non-axial radial nodes of the mode and the index I, referred to as the 
topological winding number, describes the helical structure of the wave front around a wave 
front singularity or dislocation. State vectors which are represented by a superposition of 
LG modes correspond to photons in a superposition state, with the weights of the quantum 
superposition dictated by the contribution of the Z-th angular harmonics. When the pump 
beam is a LGp° mode, under conditions of collinear phase-matching, the two-photon state at 
the output of the nonlinear crystal can be written as a coherent superposition of eigenstates 
of the OAM operator ^| that are correlated in OAM, i.e h + h = h, where l\ and I2 refer 
to the OAM eigenvalues for the signal and idler photons. A photon state described by a LG 
mode can be written as 

\lp>= J dqLG l p {q)a\q)\Q>, (1) 
where the mode function in the spatial frequency domain writes 

LG i ( OC n)-( \ 1/2 ( ^) L \i\ ( pk4) exD 



x exp <^ ilip k + i[ p - — 7r \ , (2) 



with pk and ipk being the modulus and phase, respectively, of the transverse coordinate q. 
The functions L l p (p) are the associated Laguerre polynomials and w is the beam width. 



The quantum state of the generated two-photon pair is given by 

|$ >= J dq s dqi^{q s ,qi)al{q s )al{qi)\0,0 >, (3) 

where q s ^ are the transverse coordinates in the spatial frequency domain, 1 , > is the 
vacuum state, and a\ and a\ are creation operators for the signal and idler modes. Under 
conditions of collinear propagation of the pump, signal and idler photons, the mode function 
®(<ls, Qi) can be written in the paraxial approximation [2l| as $(g s , qt) = E (q s +qi)W(q s —qi), 
where normalization of the state requires / dq s dqi\$(q s , qi)\ 2 = 1. Here E is the normalized 
field distribution of the pump beam in the spatial frequency domain, and the function W, 
which comes from the phase matching condition in the longitudinal direction, is given by 



W{qs > Qi) = ^% SmC (^4^) 6XP (^^4^) • (4) 

One can decompose the quantum state |\P > in the base of the eigenstates of the OAM 
operator, as 

l*>=EEOi.w | ^>> ( 5 ) 

h,pi h,P2 

where (h,Pi) correspond to the signal mode, (l 2 ,p 2 ) correspond to the idler mode and the 
amplitude C l ^ 2 writes 

C puP2 = J dqsdqMQs,qJ [LG l p \(q s )]* [LG 1 ^)]* . (6) 

The weights of the quantum superposition are given by Pp\fp 2 = '-J 2 , which gives the 
value of the joint detection probability for finding one photon in the signal mode (h,pi) 
and one photon in the idler mode (l 2 ,p 2 ). The two-photon state can be also characterized 
by the amplitude A(xi, z\, t\\ x 2 , z 2 , t 2 ) =< 0, 0\EiE 2 \^/ >> a ^ positions (xi, z%) and (x 2 , z 2 ) 
and times t\ and t 2 , where E + (x,z,t) is the electric field operator. If one makes use of 
he par axial approximation for the electric field operator of a photon travelling in vacuum 
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211 ]. the amplitude can be written as 
A(x 1 ,z 1} t 1 ;x 2 ,z 2 ,t 2 ) = J2J2 C^ 2 LGl\(x 1 ,z 1 )LG l ^(x 2 ,z 2 )exp(-iu(t 1 +t 2 )). (7) 

h,Pi h,P2 

Such amplitude might be employed, e.g., to draw analogies between the spatial structure of 
the two-photon quantum states and the properties of the corresponding incoherent classical 
radiation (see jia. I^). 



The amplitudes C p ^ 2 depends on two normalised parameters: the normalised pump beam 



width w p = Wp/^JXpL and the normalised beam width of the LG modes, Wq = wo/^\ p L, 
where X p is the wavelength of the pump beam in vacuum. The assumption that the diffrac- 
tion length (Ld) of the pump beam is much smaller than the length of the crystal (L « La), 
requires the condition im p w p » 1 to be fulfilled, with n p being the refractive index of the 
pump beam inside the crystal. Here we always consider w p > 1. For a typical value of the 
pump wavelength X p = 0.4/im, and a crystal length of L — 1 mm, a pump beam width of 
w p ~ 20/iin correspond to a normalized value of w p ~ 1. 

In Figure 1 we plot the contribution to the mode decomposition of the quantum state |\& > 
of all the harmonics with the same value of the indexes h and i.e. Pi lt i 2 = J2p^, P2 =o ^piiV 
In all cases cases the pump beam is a Gaussian mode. The distribution Pi lt i 2 depends only 
on the normalised width w p . Thus, for a given material, the OAM distribution of the state 



1^ > depends on the ratio w p /^X p L, a result which gives an important scaling rule. In 
Figure 2(a) we show the dependence of Pi ly i 2 on w p for different values of l\ and I2, with 
h + h = 0. One observes that the spiral bandwidth increases with w p , therefore it can be 
made larger by increasing the pump beam width or by decreasing the crystal length. Notice 
that by doing so the coupling efficiency in the detection state is also modified (see, e.g., |22j|). 
Although the total contribution Pi 1% i 2 of modes with a given index l\ and I2 depends only on 
w p , the actual amplitude distribution C^'** depends also on w . In Figure 2(b) we show the 
contributions of the modes with different index p\ and P2 for w p = 1 and ^1 = ^2 = 0, which 
is representative of a general case. We plot Cp^ ax = Sp™f>2=o Pp^p\ as a function of p ma x- 
The plot reveals that there is an optimal value of wo for which the contribution of Pq ' is 
maximum, and only a few modes with different index p make significant contributions to 
the total weight. Away from the optimal wq value, a large number of modes are required to 
represent the quantum state. 

In most applications that make use of the OAM of the photons, one projects into a 
subspace of the complete Hilbert space that describes the mode function of the photon (e.g., 
[lfij]). This implies considering only a fraction of the mode space. Let us consider only modes 
with pi = p 2 = 0, and thus thereafter we define \h,h >= \h,Pi = 0;/ 2 ,p 2 = >. In this 
subspace, the two-photon state can be written as 

00 00 

l*>= E E c l f\k,i 2 >. (8) 

ii =— 00 l 2 =~oo 
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Since the amplitude C x q 2 depends on both normalised parameters w p and w , the expansion 
given by Eq. (JHJ) depends on the base of LG modes considered. In general, calculation of 
the amplitudes C^'^ requires a 4-dimensional integration in the spatial frequency domain. 
However, for a pump beam with a spatial field profile corresponding to a LGq mode, we 
were able to obtain the value of the amplitude C 1 qq 2 is analytic form. Namely, 

<*f = ^ r ( H+m+m + /) ^ i (9) 
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xr(n) sin [ ntan -1 — 57— 1 (10) 

for h ■ h > 0. In these expressions, A stands for 

A = (-l)N-W-N ( VM f^V ( 7° ) \ (11) 

In both Eqs. Q and (jl(J|) . the relation l Q = li + Z 2 holds. It is worth stressing that Eqs. 
(jnj)- (|TTj) provide the analytical expression of the quantum state under general conditions in 
terms of the pump beam and the crystal properties. Moreover, this expansion corresponds to 
performing a Schmidt decomposition for continuos variables in this subspace, thus allowing 



to calculate the entropy of entanglement. 
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In Figure 3, we plot the expansion in OAM eigenstates of the two-photon state for pump 
beams with l = 0, 1,2 and w p = w = 1. In the three cases represented, the subspace we 
are considering represents more than 40% of the corresponding complete Hilbert space. The 
spiral bandwidth of the expansion in terms of states > increases with l . Equations 

and (fTUj) show that, for a given pump beam, the phase of the amplitude Cqq 2 changes 
for each mode with l\ ■ I2 < 0, but it is the same for all modes with li ■ I2 > 0. 

The concept can be extended to the general case of a pump beam whose spatial field 
distribution is a coherent superposition of LG modes, so that E (x) = J2 m C m LG™(x), 
with m = 0,2,4..M. The LG modes LG™(x) are the corresponding expression in the 
spatial domain of the LG modes given by Eq. (J2J). The energy flow / of the pump beam 
can be written as [15] / = 2e cn p Y^m \C m \ 2 - If we restrict to the subspace spanned by 



the states of the form \l,l >, one generates an entangled two-photon state of the form 
|\& >~ J2n=o ln\n,n >. By making use of Eqs. © and (fTUj) . one obtains that the amplitudes 
7„ of the quantum state generated when the crystal is pumped by the superposition of LG 
modes are 

2w„ 



In 
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C 2n , (12) 



This expression reveals that properly tailoring the spatial features of the beam that pumps 
the down-converting crystal allows generation of optimal quantum states for different quan- 
tum information protocols, such as the entanglement enhancement required for obtaining 
maximally entangled states. j^J One approach for the realization of entanglement concen- 
tration is to perform mode filtering operations on the two-photon state that is generated at 
the output of the crystal j^l- m contrast, Eq. (fT2j) shows that the required quantum state 
could be obtained through appropriate engineering of the spatial properties of the pump 
beam 2(|. It also shows that the phase of the state \n,n > corresponds to the phase of 
the corresponding mode of the classical pump beam, arg{C2 n }, a result consistent with the 
experimental observation by Ou et. aL that the down-converted photons carry information 
about the phase of the pump beam |27j ]. 

In conclusion, we have obtained and analyzed the detailed quantum spatial structure 
of the two-photon entangled states generated in parametric down-conversion in terms of 
the eigenstates of the orbital angular momentum operator. We put forward the related 
concept of quantum spiral bandwidth, and showed its dependence on the pump beam and 
down-converting crystal. This allows, e.g., to define an effective finite Hilbert space where 
entanglement takes place. [3] Engineering such quantum spatial bandwidth should be a key 
tool to optimize multidimensional quantum information schemes, for example to increase 
the resolution of two-photon imaging using entangled photons (e.g., J^]) and to enhance 
the efficiency of relevant multidimensional quantum cryptography protocols (e.g., js^|). 

This work was supported by the Generalitat de Catalunya and by the Spanish Govern- 
ment through grant BFM2002-2861. 
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Figure captions 

Fig. 1: Mode distributions Pi 1: i 2 for different values of the parameter w p . (a) w p = 1, (b) 
w p = 2.5 and (c) w p = 5. The pump beam is a gaussian mode {l p = 0), so that l± + 1 2 — 0. 
The represents the value of l\. 

Fig. 2: (a) Value of the weight P Jljj2 for (Z x , Z 2 ) = (0,0), (1, —1) and (2, —2), as a function 
of the normalized pump beam width w p . The labels show the value of (hjh)- The pump 
beam is a gaussian mode, (b) Weight C®'® ax as a function of the index p rnax for different 
values of w , as shown in the labels. In all cases w p — 1. 

Fig. 3: Mode decomposition in the subspace \h,pi = 0; Z 2 ,P2 = > for several pump beam. 
In (a), (c) and (e) we plot the weight Pq 1 ^ 2 of the mode distribution , and in (b), (d) and 
(f) we plot the phase arg{C l Q Q 2 } of each mode, (a), (b): l = 0; (c), (d): l — 1 and (e), (f): 
Z = 2. 
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